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Abstract 

Understanding the interaction between sexual and natural selection within variable 
environments is crucial to our understanding of evolutionary processes. The hand- 
icap principle predicts females will prefer males with exaggerated traits provided 
those traits are indicators of male quality to ensure direct or indirect female bene- 
fits. Spatial variability in ecological factors is expected to alter the balance between 
sexual and natural selection that defines the evolution of such traits. Male and 
female blackspotted topminnows (Fundulidae: Fundulus olivaceus) display promi- 
nent black dorsolateral spots that are variable in number across its broad range. We 
investigated variability in spot phenotypes at 117 sites across 13 river systems and 
asked if the trait was sexually dimorphic and positively correlated with measures of 
fitness (condition and gonadosomatic index [GSI] ). Laboratory and mesocosm ex- 
periments assessed female mate choice and predation pressure on spot phenotypes. 
Environmental and community data collected at sampling locations were used to 
assess predictive models of spot density at the individual, site, and river system level. 
Greater number of spots was positively correlated with measures of fitness in males. 
Males with more spots were preferred by females and suffered greater mortality 
due to predation. Water clarity (turbidity) was the best predictor of spot density on 
the drainage scale, indicating that sexual and natural selection for the trait may be 
mediated by local light environments. 



Introduction 

Sexual selection is known to be a potent evolutionary force. 
Understanding how various sexual selection mechanisms 
work within the context of variable environments is crucial 
to our understanding of evolutionary processes (Candolin 
and Heuschele 2008; Walsh and Reznick 2009; Ingleby et al. 
2010). The handicap principle predicts that females will pre- 
fer males with exaggerated traits that are indicators of genetic 
quality (Iwasa et al. 1991; Iwasa and Pomiankowski 1999; 
Cotton et al. 2009). Male-exaggerated traits must be indi- 
cators of increased heritable fitness if females are to benefit 
indirectly through increased offspring fitness. The cost of ex- 
aggerated male traits may be extensive with natural selection 
mediating the trade-off between male fitness and expression 
of exaggerated traits (Rosenthal et al. 200 1 ) . Spatial and tem- 
poral environmental variability can alter this balance between 



sexual and natural selection, preventing populations from be- 
coming fixed by persistent directional selection (Moller and 
Alatalo 1999; Sharma et al. 2010; Vergara et al. 2011). Local 
environmental conditions can dramatically affect the efficacy 
of sexual signals (Fuller and Travis 2004; Millar et al. 2006; 
Vergara et al. 201 1) and survival rates of males expressing ex- 
aggerated traits (Houde and Endler 1990; Millar et al. 2006; 
Gordon et al. 2011). Empirical data on broadly distributed 
sexually dimorphic species found in a range of ecological 
conditions is needed to better understand the dynamics of 
natural and sexual selection within a meaningful ecological 
context. 

Freshwater fish have served as model study systems for a 
number of fundamental questions in evolutionary ecology. A 
major advantage of these systems is that streams feature linear 
and predictable gradients in ecologically important bio tic and 
abiotic factors (Vannote et al. 1980). In many cases, natural 
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barriers within or among river systems form replicate evolu- 
tionary units. In Trinidad killifish (Rivulus hartii), predation 
pressure and resource availability correlated with stream gra- 
dients shape the evolution of a number of life-history traits 
(Fraser et al. 1995; Reznick et al. 2001; Walsh and Reznick 
2009, 2010). Predation pressure and flow regimes influence 
the evolution of body shape and a number of correlates of 
performance in Gambusia affinis (Langerhans et al. 2004; 
Langerhans 2008; Langerhans and Reznick 2010). Studies 
of sexual and natural selection within an ecological context 
have shown female choice, light environments, and preda- 
tion pressure to interact in driving the evolution of male 
phenotypes (Houde and Endler 1990; Endler 1993; Millar 
et al. 2006). Millar et al. (2006) linked the relative number 
of different color spots to water clarity and the physiological 
constraints of predator vision systems in contrasting envi- 
ronments. While these and other well-studied systems serve 
as foundations to our understanding of the evolutionary pro- 
cess, they are often limited in scope. Predation regimes are 
typically a dichotomous presence or absence and the distribu- 
tion of the species studied are either limited or only studied 
in detail in a limited area (Krebs and Bell 2011). We com- 
bined laboratory experiments of mate choice and predation 
pressure with fieldwork conducted along stream gradients 
across 13 drainages to address questions regarding observed 
variability in the expression of a sexually dimorphic trait. 

The genus Fundulus contains 38 extant species that in- 
habit a variety of habitats in North and Central America 
(Fuller et al. 2007; Whitehead 2010). The majority of species 
in the group are sexually dimorphic with individual species 
expressing a variety of sexually dimorphic traits such as ver- 
tical bars, lateral stripes, color patterns on the body or fins, 
and various color spotting patterns. Within the group, ge- 
netic and plastic (linked with water clarity and transmission 
of ultraviolet [UV] ) effects influence vision systems and color 
phenotypes of Lucania goodei (Fuller and Travis 2004; Fuller 
et al. 2005). The group selected for this study is the widely 
distributed Fundulus notatus species complex, which occurs 
throughout most of the central United States. The F. notatus 
species complex is comprised of the Blackstripe Topmin- 
now (F. notatus), the Blackspotted Topminnow (F. olivaceus), 
and the Broadstripe Topminnow (F. euryzonus) that all share 
similar niches (Thomerson and Woolridge 1970; Blanchard 
1996) and hybridize in contact zones found throughout 
their ranges (Thomerson 1967; Setzer 1970; Duvernell et al. 
2007; Schaefer et al. 2011). Throughout most of their range, 
F. olivaceus tends to occupy high gradient, clear streams, while 
F. notatus dominates more turbid backwater habitats of larger 
rivers and prairie streams (Howell and Black 1981). All three 
species in the F. notatus complex have a distinct lateral stripe 
that extends the length of the body along with black spots 
typically above the stripe. The number and intensity of these 
spots is cited as a diagnostic character for identification of 
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F. olivaceus (Lee et al. 1980; Suttkus and Cashner 1981; Ross 
2000). While spots of various colors are seen in a number 
of Fundulus, black spots are unique to the P. notatus com- 
plex. Thomerson (1966) noted variability in spot phenotypes 
across the range and hypothesized it may be a sexually di- 
morphic trait. In our own work with these species (Duver- 
nell et al. 2007; Schaefer et al. 2009, 2011), we have noted 
substantial variation in F. olivaceus spot phenotypes among 
river systems and sexes and recognized that this trait is not 
always reliably diagnostic for species identification. 

Our first objective was to quantify variability in the 
spot phenotype between males and females of the two 
more widespread species in the complex that display spots 
(F. notatus and F. olivaceus). We ask if the spot phenotype 
is sexually dimorphic and correlated with measures of con- 
dition and fitness for either species. Observed variability in 
the spot phenotype and results of the above analyses moti- 
vated us to ask detailed questions about the spot phenotype 
in F. olivaceus. Specifically, (1) do females select male F. oli- 
vaceus based on the spot phenotype? (2) is predation risk 
related to the spot phenotype? (3) what ecological variables 
are the best predictors of variability in the spot phenotype of 
F. olivaceus across its broad range? We assessed the efficacy 
of spot phenotype models at the individual, population, and 
river system level. Models included measures of key ecolog- 
ical traits hypothesized to be related to natural and sexual 
selection dynamics (predation pressure, habitat quality, and 
light environment). 

Methods 
Fish collection 

Within each of the 13 river systems we identified between 
8 and 12 sites (Table 1) for sampling. Based on museum 
records and earlier work with these species (Duvernell et al. 
2007; Schaefer et al. 2011), sampling in each drainage was 
centered around the confluence of smaller tributaries (typ- 
ical F olivaceus habitat) and larger rivers (F. notatus habi- 
tat) with roughly half the sampling sites upstream in trib- 
utaries and downstream in larger rivers (Fig. 2). This en- 
sured that both species and potential areas of coexistence 
were sampled along a gradient of stream sizes. At each site, 
we collected fish using a dip net and a 6.1 x 1.2 m seine 
with 0.3-cm mesh. Dip nets were used to specifically target 
Fundulus sp. that are easily visible at the surface and effi- 
ciently netted. Fish assemblages (including a smaller number 
of both species of Fundulus) were sampled by seining all 
available habitats over a 100-150 m reach. For each indi- 
vidual Fundulus sp. collected, a tissue sample was preserved 
for later genetic identification (Schaefer et al. 2011) before 
the fish was placed in an individually labeled 50-mL coni- 
cal tube containing 10% formalin. Fish assemblage samples 
were preserved in 10% formalin and later transferred to 70% 
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Table 1. Number of sites sampled, number of adult (>32 mm SL) Fundulus analyzed, and a summary of the range of ecological variables measured 
at sites within each river system. Predation pressure is expressed as the log abundance of piscivorous fish and size represents cumulative drainage 
area. Names of the river systems as in Figure 1 . 



River system 


Sites 


F. notatus 


F. olivaceus 


Predation (log abundance) 


Turbidity (NTU) 


Canopy (%) 


Size (km 2 ) 


Amite 


10 


0 


167 


0.0-1.2 


4.6-20.1 


0.0-80.0 


7-318 


Black 


9 


9 


202 


0.3-1.1 


1.7-52.4 


1.7-52.4 


32-21,397 


Elk 


7 


35 


82 


0.5-0.8 


7.2-20.7 


10.0-80.0 


61-3,97 


Illinois 


7 


0 


153 


0.3-1.7 


4.11-9.73 


1.7-38.3 


638-3840 


Kiamichi 


8 


0 


235 


0.3-1 .0 


3.1-27.6 


3.3-25.0 


12-1806 


Little 


8 


36 


111 


0.3-0.9 


2.44-12.2 


1.7-18.3 


486-2996 


Neches 


8 


53 


117 


0.0-0.6 


9.4-24.2 


0.0-73.3 


45-20,500 


Saline 


12 


102 


160 


0.0-1.1 


1.5-189.0 


0.0-80.0 


4-498 


Pascagoula 


9 


86 


75 


0.0-1.1 


6.3-15.8 


3.3-60.0 


1738-21,177 


Pearl 


10 


80 


159 


0.0-0.5 


3.8-42.5 


0.0-60.0 


35-16,555 


Sabine 


9 


5 


101 


0.3-0.8 


6.2-123.0 


0.0-48.3 


45-21,530 


Spring 


10 


42 


165 


0.0-0.8 


5.9-26.1 


3.0-35.0 


337-6318 


Tombigbee 


10 


72 


113 


0.6-1.2 


3.8-14.0 


1.7-70.0 


31-6052 


Totals 


117 


520 


1840 


0.0-1.7 


1.5-189.0 


0.0-80.0 


4-21,530 




Figure 1. Photographs of large adult male Fundulus olivaceus ([A] from 
the Little River and [B] from the Saline River drainage) and F. notatus ([C 
and D] both from the Pearl River drainage) demonstrating the range of 
spot phenotypes seen in the species. 

ethanol, identified to species, enumerated, and deposited into 
The University of Southern Mississippi Ichthyological Col- 
lection (http://ichthyology.usm.edu/usm/). At each site, we 
measured turbidity (NTU, HACH 2 1 00 turbidity meter) and 



estimated canopy cover (mean of three estimates taken along 
transects at the upper, middle, and lower portions of sampled 
reach of stream). 



Spot density and individual condition 

For each Fundulus sp., we counted the total number of visible 
dorsolateral spots on the left side of the body and measured 
standard length (SL), wet weight (mass), gonad mass, and 
body mass without gonads or gastrointestinal track (evis- 
cerated mass). For analyses, we used only adults (>32 mm 
SL) for which sex could be determined reliably. Because the 
number of spots was positively correlated with SL (r 2 = 0.43, 
P < 0.0001, df = 1,2573; spot count = 1.71 x SL - 21.5), 
we controlled for body size by dividing the number of spots 
by SL (hereafter spot density). To test and control for allom- 
etry, we first tested for it by regressing spot density and SL 
in each species-sex combination and then included SL as a 
covariate in analyses where appropriate. Body condition was 
assessed using Fulton's Condition Index (1000 x eviscerated 
mass/SL 3 ) (Sutton et al. 2000) and reproductive condition as- 
sessed as gonadosomatic index (GSI, ratio of gonad mass to 
eviscerated mass). Individuals were assigned to species using 
diagnostic nuclear intron markers at five or more indepen- 
dent loci following Schaefer et al. (2011). Any individuals 
who were not homozygous at all five diagnostic loci were 
excluded from analyses. We used a two-way analysis of co- 
variance (ANCOVA) (SL as covariate) to test for differences 
in spot density among species and sex. For each species and 
sex combination, we correlated spot density (SL as a covari- 
ate) with body condition and spot density with GSI to assess 
whether spot density was an indicator of condition or fitness. 
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Spot density and male reproductive success 

To assess the influence of spot density on reproductive suc- 
cess, we measured the proportion of offspring sired in labo- 
ratory mate choice trials. Male and female F. olivaceus were 
collected from creeks within the Pascagoula River drainage 
and housed in the University of Southern Mississippi wetlab 
facility. Males were initially screened at a single polymor- 
phic microsatellite locus (FhATG-B103) (Adams et al. 2005; 
Duvernell et al. 2007; Gutierrez 2010). For each trial, four 
males with varying spot phenotypes whose offspring could 
be unambiguously identified with the single locus were se- 
lected, anesthetized, and photographed to count spots and 
measure SL. We calculated spot density as described above 
and ranked the four males in spot density. The four males and 
either one or two females were placed in 200-L round pools 
with gravel substrate and four acrylic mops in various loca- 
tions of the pool that served as spawning media. Spawning 
mops were checked daily and a trial was considered complete 
when viable clutches were collected on four consecutive days. 
Eggs were incubated in smaller plastic cups for 10-20 days 
(Vigueira et al. 2008) and hatched larvae were preserved for 
genotyping after hatching. Trials produced a minimum of 
10 eggs and when a larger number of eggs (80 or more in 
three trials) were produced roughly 25% of hatched larvae 
were selected at random for genotyping. After trials, males 
were again anesthetized and photographed to ensure spotting 
patterns did not vary significantly over the length of the trial. 
There was minimal variability in spot patterns and fish could 



easily be identified when comparing before and after photos. 
For a total of 13 trials (seven with one female and six with 
two females), we quantified the proportion of larvae sired by 
each of the four males. We used ANCOVA to test for differ- 
ences (across all trials) in logit-transformed proportions of 
offspring sired among the four rank spot densities with SL as 
a covariate. 

Spot density and predation 

To assess the influence of spot density on predation pres- 
sure, we measured survivorship when exposed to a preda- 
tor (spotted bass, Micropterus punctulatus) (Knight and 
Gido 2005). Predation trials were conducted in a series of 
24 outdoor mesocosms designed to mimic small stream habi- 
tat (Matthews et al. 2006). Mesocosms consisted of a pool 
(1.83 m diameter, 0.6-m deep with mixed sand/gravel sub- 
strate and woody debris as structure) and a riffle (0.9 in- 
king, 0.3-m wide, and 15-cm deep with gravel/cobble sub- 
strate) with recirculating pumps and a groundwater source. 
Fundulus olivaceus and M. punctulatus (135-210 mm SL) 
were collected from local tributaries in the Pascagoula River 
drainage and housed in three mesocosms (one for M. punc- 
tulatus and for each F. olivaceus sex) until trials began. For 
each trial, eight F. olivaceus (male and female trials done sepa- 
rately) were anesthetized and photographed to calculate spot 
density as described above. Because prey size is expected to 
strongly influence predator behavior, trials consisted of fish 
of similar sizes (mean SL of 54.5 mm for all trials, mean size 
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range within trials was 4.5 mm with a maximum in one trial of 
11 mm). Fish were introduced to mesocosms and allowed to 
acclimate for at least 30 min before a predator was introduced. 
Mesocosms were checked every 24 h to assess how many fish 
remained and trials ended when roughly half of the prey had 
been consumed (typically two to six days). Trials in which all 
fish were consumed within the first 24 h were excluded from 
analyses. The surviving fish were photographed and photos 
matched to pretrial photos. Six control trials were run with 
no predators to ensure that all fish were recovered and all 
mortality in experimental trials was due to predation. The 
eight fish in each trial were categorized as high or low spot 
density (four highest and lowest spot densities, respectively) 
and we tested for differences in logit-transformed survivor- 
ship between high and low spot density trials with paired 
f-tests for male and female trials. 

Individual, population, and drainage level 
variability in spot density 

We used Akaike's Information Criteria (AIC C ) for small sam- 
ples size to compare the predictive power of candidate models 
for male and female F. olivaceus spot density at the individual, 
site, and drainage level (Anderson et al. 2000 ) . Predictive vari- 
ables (Table 2) included predation pressure (log abundance of 
all piscivorous species collected [Lepisosteus sp., Micropterus 
sp., Esox sp., Lepomis cyanellus, and L. gulosus]), water clar- 
ity (turbidity), reproductive investment (GSI), condition 
(Fulton's condition index), canopy cover, and drainage (river 
system). For individual male and female models, condition 



and GSI represented measures of individual fitness. For pop- 
ulation level models, GSI and condition were averaged for 
each sex at each site and served as a measure of habitat quality. 
For river system models of spot density, these variables along 
with predation pressure, canopy cover, and turbidity were av- 
eraged across all sites within each river system. The drainage 
variable was not included in river system models. Candidate 
models included a null (no variables), each one of these vari- 
ables alone and all possible combinations of two-way models. 
Individual, site or drainage mean SL was included in all mod- 
els except for the null to control for allometry. Three way and 
higher interactions were not included in analyses due to diffi- 
culties with interpretation. Models with low dAIQ and high 
Akaike weights ( w;) have the best combination of parsimony 
(few parameters) and predictive power. We only interpreted 
models with dAIQ <2.0 and w, score > 10% of the high- 
est Wi score as meaningful (Anderson and Burnham 2002). 
All analyses were performed in the R statistical language 
(R Development Core Team 2010). 

Results 

Variability in spot density among sex and 
species 

We collected data on 2575 adult Fundulus from 117 sites in 
the 13 river systems. The number of spots ranged from 0 (sev- 
eral male and female F. notatus from the Pascagoula River) to 
over 150 (male P. olivaceus from tributaries in the Spring 
and Black River systems) (Fig. 1). Spot densities showed 
positive allometry for male F. olivaceus (SL-spot density 



Table 2. Variables used in AIC modeling of individual, population, and river system level models of male and female spot density in 
Fundulus olivaceus. For individual models, individual reproductive investment (GSI) and condition (Fulton's condition index) values were used. For 
population and river system level models, those variables were averaged as indicators of habitat quality. Other site-specific variables (predation pres- 
sure, water clarity, and canopy cover) were averaged for river system models. At each of the three levels, candidate models included all variables alone 
and all possible combinations of two-variable models with interaction terms. Standard length (individual, site, or drainage mean) was a covariate in all 
models to control for allometry. 



Variable 



Measure 



Hypotheses 



Predation pressure 

Reproductive investment or habitat 
quality 



Condition or habitat quality 

Water clarity 
Canopy cover 
Drainage 



Log abundance of piscivores in local Increased predation pressure selects against increased spot density, 
fish community 

Spot density is an indicator of fitness on the individual level. On the 
site (population) and drainage level, more favorable habitat will 
have healthier fish that have larger gonads and display more 
spots. 

Spot density is an indicator of condition on the individual level. On 
the site (population) and drainage level, more favorable habitat 
will have healthier fish that are in better condition and display 
more spots. 

Decreased water clarity reduces sexual selection for increased spot 
density. 

Decreased light availability reduces sexual selection for increased 
spot density. 

Spot density is not under strong selection and patterns will best be 
explained by random differences among drainages (drift). 



Gonadosomatic index 

Fulton's condition index 

Turbidity (NTU) 

Canopy cover (%) 

13 river systems sampled 
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F. notatus 

n.270m,250l 



F. olivaceus 

n-89l m. 949 I 



Figure 3. Mean (±2 SE) spot density for male and female Fundulus 
notatus and F. olivaceus. 



slope = 0.021, r 2 = 0.213, P < 0.001) and negative allome- 
try for female F. olivaceus (slope = -0.003, r 2 = 0.008, P < 
0.001) and male F. notatus (slope = -0.012, r 2 = 0.046, both 
P < 0.001), but no pattern of allometry in female F. notatus 
(slope = -0.005, r 2 = 0.009, both P < 0.146). There was a 
significant difference in spot density among the two species 
(F. olivaceus mean = 1.34 ± 0.011 SE;P. notatus mean = 0.798 
± 0.018 SE; J 7 !, 2355 = 74.3, P < 0.001, r)\ = 0.22) and an in- 
teraction (Fi i23 55 = 22 3.7, P < 0.001, r)\ = 0.07) between 
species and sex (Fig. 3). This interaction was due to male 
F. olivaceus having significantly higher spot density than fe- 
males (-F U 837 = 122.6, P < 0.001, r]\ = 0.31) while the sexes 



did not differ in F. notatus (F^sn 
0.01). 



0.07, P < 0.79, 



Spot density and individual condition 

Spot density was positively correlated with condition in male 
F. olivaceus (F lM9 = 23.7, P < 0.001, r 2 = 0.23) and male 
F. notatus (F lj26 8 = 24.62, P < 0.001, r 2 = 0.13). Spot density 
was positively correlated with condition in female F. olivaceus 
{F lM7 = 31.69, P < 0.001, r 2 = 0.04) and female F. notatus 
(F h24g = 20.7, P < 0.001, r 2 = 0.12) (Fig. 4). Spot den- 
sity was positively correlated with GSI in male P. olivaceus 
(Pi, 889 = 96.70, P < 0.001, r 2 = 0.27) and male F. notatus 
(F 1,268 = 20.70, P < 0.001, r 2 = 0.12) (Fig. 5). Spot den- 
sity was negatively correlated with GSI in female F. olivaceus 
(i? 1947 = 8.90, P < 0.001, r 2 = 0.02) but positive in female 
F. notatus (Fi j2 as = 16.90, P < 0.001, r 2 = 0.08). 

Reproductive success and predation pressure 

Within mate choice trials, the range (difference between first 
and fourth ranked individuals) of spot densities for males av- 
eraged 0.68 (minimum of 0.31, maximum of 1.10). There 
was a significant difference in male reproductive success 
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Figure 4. Relationship between spot density and Fulton's condition in- 
dex in female (top panel) and male (bottom panel) Fundulus notatus and 
F. olivaceus. Lines represent best fit linear regressions. 

among the four spot density ranks (^3,47 = 3.39, P < 0.025, 
rj 2 , = 0.16). Males that ranked first in spot density sired on 
average 53.8% of offspring that was significantly greater than 
males that ranked second (15.4%, Tukey multiple compari- 
son P < 0.021) or third (7.6%, P < 0.006), but not different 
than males that ranked fourth (23.1%, P < 0.063). There 
were no pairwise significant differences in male reproduc- 
tive success among second, third, or fourth ranked males 
(Fig. 6A). Within predation trials, the range (difference be- 
tween first and fourth ranked individuals) of spot densities 
averaged 0.82 (minimum of 0.31, maximum of 1.31). The 
rate of survival in predation trials was significantly lower for 
high spot density males compared to low spot density males 
(43.8 vs. 65.6% survival, t = -2.73, P < 0.016). There was no 
difference in survivorship between high and low spot density 
females (50.0 vs. 59.4%, t = -0.38, P < 0.72) (Fig. 6B). 

Models of F. olivaceus individual, 
population, and river system spot density 

Individual spot density was best predicted by a single model 
in male (Drainage x GSI) and female (Drainage x Canopy; 
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Figure 5. Relationship between spot density and gonadosomatic index 
in female (top panel) and male (bottom panel) Fundulus notatus and 
F. olivaceus. Lines represent best fit linear regressions. 



Table 3) F. olivaceus. Both models contained drainage as 
a variable and predicted between 39% and 53% of the 
variability among individuals. On a population level, 
drainage was again a component of both interpretable mod- 
els. Average male spot density at sites was best predicted by 
a single model (Drainage) and average female spot density 
was best predicted by the Drainage x GSI model. Differences 
among drainages were substantial with nearly 63% of the 
variation at the site level explained by drainage differences 
(Table 3; Fig. 7). On the river system scale, both male and 
female mean spot density was best predicted by Turbidity 
alone and for females the Turbidity x GSI model was also 
interpretable. Across drainages, Turbidity alone was a better 
predictor of male than female spot density (Table 3), but the 
slopes of the relationship for males and females were very 
similar (Fig. 8). 

Discussion 

Patterns of variability in the spot phenotype for F. olivaceus 
were consistent with the handicap principle. In F. olivaceus, 
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Figure 6. Proportion of offspring sired by rank male spot density (top 
panel) and mean survival rate for males and females with high or low 
spot densities. 



spot density was an indicator of male fitness and females pre- 
ferred males with more spots in mate choice trials involving 
fish from the Pascagoula River drainage. Predation pressure 
as a putative cost to higher spot density was supported in 
mesocosm predation trials (fish from the Pascagoula River 
drainage) but not in field data as predator abundance was not 
part of interpretable models at any level. In many well-studied 
systems where predation pressure is consistently shown to 
influence male phenotypes (Rosenthal et al. 2001; Langer- 
hans et al. 2004; Millar et al. 2006), predation pressure is 
effectively depicted as a dichotomous presence or absence 
(Krebs and Bell 20 1 1 ). In this system, predation is more com- 
plex as predators (large-bodied piscivores) were in all likeli- 
hood ubiquitous. The sampling gear used (seine) has been 
shown to preferentially select small-bodied species (Gido 
et al. 2009), meaning we likely undersampled predators, espe- 
cially in larger river habitats. The expectation ( Vannote et al. 
1980) is that the abundance and size of predators will increase 
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Table 3. Model complexity (K), AIC C scores, criteria for model selection (dAIC c , w), and adjusted r 2 for individual, population, and river system level 
models of male and female spot density in Fundulus olivaceus. Only interpretable models (dAIC c <2.0 and iv ; score >10%) are presented. Standard 
length was a covariate in all models to control for allometry. 



Model 



K 



AIC C 



dAICc 



W, 



Individual 
Male F. olivaceus 
Drainage x GSI 
Female F. olivaceus 
Drainage x Canopy 
Population 
Male F. olivaceus site means 
Drainage 
Female F. olivaceus site means 
Drainage x GSI 
River system 
Male F. olivaceus 
Turbidity 

Female F. olivaceus 
Turbidity 
Turbidity x GSI 



28 
28 

15 

27 



455.2 
230.1 

7.5 
725.0 

-3.3 

101.0 
101.8 



0.0 
0.0 

0.0 
0.0 

0.0 

0.0 
0.8 



0.99 
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0.987 
0.497 

0.892 

0.549 
0.377 
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0.392 

0.626 
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0.771 
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Figure 7. Mean spot density (±1 SE) by river system for mature male 
Fundulus notatus and F. olivaceus. 
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Figure 8. Relationship between mean turbidity and spot density among 
the 13 river systems for adult male and female Fundulus olivaceus. Values 
are the means among all sites in each system ±1 SE. 



with increasing stream size. Turbidity is expected to (and did) 
follow the same pattern of increasing with stream size within 
a river system. Thus, it is difficult in some respects to dis- 
entangle the potential influence of turbidity and predation 
pressure within river systems in this study. However, there 
were clear differences in turbidity among river systems that 
correlated strongly with river system level differences in spot 
density. There is no reason to expect (from our data or what 
is known of North American fish community composition) 
significant drainage level differences in predator abundance 
(Matthews 1998). 

Assuming predation pressure is generally consistent across 
river systems and the abundance of large-bodied piscivores 
increases in downstream habitat (Vannote et al. 1980), there 



are a number of indicators that predation has played an im- 
portant role in the evolution of F. notatus and F. olivaceus. The 
expectation (Reznick and Endler 1982; Reznick et al. 2001; 
Krebs and Bell 20 1 1 ) is that predation pressure will select for a 
suite of life-history traits that are consistent with observed dif- 
ferences between these two species. Fundulus notatus (higher 
predation habitats) is smaller (Lee et al. 1980), reproduces 
earlier at a smaller size (J. F. Schaefer, unpubl. data), pro- 
duces a larger number of smaller offspring (Vigueira et al. 
2008) and invests more energy in reproduction (Fig. 4). It 
is also intriguing that within Fundulus, the black spot phe- 
notype is unique to these species and evolved as a sexually 
dimorphic trait in upstream habitats that feature clearer wa- 
ter and presumably lower predation pressure. 
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Genetic versus plastic contributions 

We do not have data to directly address the relative contri- 
bution of genetic and environmental influences on the spot 
phenotype. While both mechanisms are likely involved with 
expression of this trait (Fuller and Travis 2004), a number of 
observations would be consistent with a fairly limited adult 
plasticity. First, spot phenotypes did not change over short 
periods of time (4-21 days) in predation (conducted in an 
outdoor setting with natural food, light, and photoperiod) 
or mate choice (conducted indoors with artificial food and 
lighting) experiments. In all cases, spot patterns acted like 
a fingerprint and individual fish could be unambiguously 
matched through time. For predation trials, we had two au- 
thors independently match individuals using the same digital 
images. Second, spot phenotypes (unlike fin color) are gener- 
ally not seasonally variable. While not quantified, our expe- 
rience working with these species is that individuals caught 
in early spring or late fall (outside the breeding season) still 
express the spot phenotype. This general consistency in adult 
spot patterns is likely how this trait came to be regarded as 
diagnostic (Suttkus and Cashner 1981) among the species. 
While this indicates that adult phenotypes may be somewhat 
fixed, the trait clearly changes allometrically (see above) and 
there may be a considerable developmental plasticity ( West- 
Eberhard 2003) component. A heritability study is planned to 
directly address the genetic contribution to this phenotype. 

Natural and sexual selection 

The exact mechanisms that lead to either increased preda- 
tion or reproductive success of higher spot density fish is not 
known. Species in the F. notatus species complex are counter 
shaded (unpigmented below and darker above) and the broad 
lateral stripe extending from the tail through the eye (unique 
among Fundulus) is most likely an example of disruptive col- 
oration (Burtt 1981; Stevens and Merilaita 2009). Disruptive 
color patterns feature areas of high contrast to break up body 
shape in which case high spot density may make males more 
conspicuous by reducing contrast or outlining the body. It is 
also possible that the spots are not the direct cause of increased 
predation risk. Larger reproductive investment may merely 
be correlated with spot density while behavioral changes in- 
crease exposure to predators and reproductive success. In 
documenting mating behaviors, substantial male-male and 
even female-female aggression was observed along with male 
lateral displays to females (Gutierrez 2010). The potential im- 
portance of male dominance was also evident in two mate 
choice trials where one of the successful males was a sub- 
stantially larger fish with relatively few spots (ranked fourth 
in spot density). We did not collect behavioral data as part 
of these trials, but our results are consistent with females 
preferring size and high spot density or male-male domi- 
nance resulting in larger males having more access to females 
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(Fig. 4). It should also be pointed out that the arenas where 
mate choice trials took place were fairly small (1-m diameter 
pools), which would exacerbate the importance of aggressive 
interactions. Regardless of the mechanisms involved, females 
preferred increased spot density that clearly had a fitness cost 
to males in the form of increased predation risk. 

Geographic patterns 

One of the more appealing aspects of this experimental system 
is the species broad distribution across ecologically disparate 
river systems (Schaefer et al. 2011). Differences among river 
systems were clearly important (Fig. 6), as the drainage vari- 
able was part of every interpretable spot density model at the 
individual and population level (Table 3). Differences among 
river systems were not random (expectation if selective pres- 
sures were weak) as turbidity was strongly correlated with 
spot phenotypes on the river system level (Table 3; Fig. 7). 
The general pattern revealed was that drainages in the Coastal 
Plain (right side of x-axis in Fig. 6, Amite, Neches, Pascagoula, 
Pearl, Sabine, and Tombigbee) have lower spot densities than 
those in the Ozark (Illinois, Black, and Spring), Ouachita 
(Little and Kiamichi), or Eastern Highlands (Elk). The lone 
river system outside of the Coastal Plain with low spot den- 
sity is in the Saline River, a tributary of the Ohio River in 
the formerly glaciated Southern Till Plain of southern Illi- 
nois. The Ozark, Ouachita, and Eastern Highland tributaries 
(F. olivaceus habitat) typically feature larger (cobble and 
gravel) and more stable substrate compared to prairie or 
coastal systems where tributaries typically have finer gravel, 
sand, silt, or clay substrate. These differences likely result in 
observed turbidity patterns across river systems. We can also 
generally rule out phylogeographic influence as the cause 
of river system differences in spot density. A population 
level phylogeny (B. R. Kreiser, unpubl. data) has revealed 
relatively little structure across the range for F. olivaceus 
(distinctly different from patterns in F. notatus). Even with 
modest grouping of coastal versus inland populations (sup- 
ported by phylogeography, B. R. Kreiser, unpubl. data), rela- 
tionships between turbidity and spot density persist within 
these two groups. 

While spot expression in male F. olivaceus fits predictions 
of the handicap principle, females also express the pheno- 
type. In many well-studied sexually dimorphic fishes, genes 
controlling dimorphic color traits are often linked to the Y- 
chromosome and not expressed in females (Lindholm and 
Breden 2002; Gordon et al. 201 1). Female expression of spots 
raises the possibility that the spot phenotype in these fish has 
a different genetic architecture or control mechanism. One 
possibility is hormonal regulation, something common in 
sexually dimorphic color patterns (Parker et al. 20 1 1 ) . Hask- 
ins et al. (1961) found that testosterone treatments led to the 
expression of non-Y-linked color genes in female guppies. 



©2012 The Authors. Published by Blackwell Publishing Ltd. 



1379 



Evolution of Sexual Dimorphism in Topminnows 



J. F. Schaefer ef al. 



Hormonal regulation of spot density in F. olivaceus would 
be consistent with the observed strong male allometry (spot 
density increases as males mature) and the negative correla- 
tion between female spot density and GSI (Fig. 5) but not 
female condition (Fig. 4). The similarity in male and female 
spot densities across drainages (Fig. 8) is intriguing as a signal 
of parallel evolution of the sexes within a species, something 
generally not well studied (Hendry et al. 2006). As a whole, 
our data are consistent with natural and sexual selection act- 
ing on males driving parallel evolution of both sexes in a trait 
that is hormonally regulated. 

Data for this study were compiled over distinct linear gra- 
dients within a broad range of disparate river systems. The re- 
sult is compelling support for some fundamental processes in 
evolutionary ecology. Future work with this system is planned 
to address the underlying genetic framework and heritabil- 
ity of the spot phenotype and other suspected sources of 
variability (such as plasticity). The exact role of spots as a 
reproductive isolating mechanism among the species has not 
been tested but one would suspect some role based on mate 
choice data presented. The potential interaction between the 
local environmental gradients, expression of traits contribut- 
ing to reproductive barriers, and the structure of contact 
zones (Schaefer et al. 2011) provides ample opportunity to 
address a number of stimulating questions. 

Acknowledgments 

Funding was provided by the National Science Foundation 
(DEB no. 0716985, DEB no. 0716964, and MRI 0923063). 

References 

Adams, S. M., M. F. Oleksiak, and D. D. Duvernell. 2005. 

Microsatellite primers for the Atlantic coastal killifish, 

Fundulus heteroclitus, with applicability to related Fundulus 

species. Mol. Ecol. Not. 5:275-277. 
Anderson, D. R., and K. R Burnham. 2002. Avoiding pitfalls 

when using information-theoretic methods. ]. Wildl. Manage. 

66:912-918. 

Anderson, D. R., K. P. Burnham, and W. L. Thompson. 2000. 

Null hypothesis testing: problems, prevalence, and an 

alternative. ]. Wildl. Manage. 64:912-923. 
Blanchard, T. A. 1996. Ovarian cycles and micro habitat use in 

two species of topminow, Fundulus olivaceus and F. euryzonus, 

from the southern United States. Environ. Biol. Fish. 47:155- 

163. 

Burtt, E. H. 1981. The adaptiveness of animal colors. Bioscience 
31:723-729. 

Candolin, U., and ]. Heuschele. 2008. Is sexual selection 

beneficial during adaptation to environmental change? Trend. 

Ecol. Evol. 23:446-452. 
Cotton, S., ). Small, R. Hashim, and A. Pomiankowski. 2009. 

Eyespan reflects reproductive quality in wild stalk-eyed flies. 

Evol. Ecol. 24:83-95. 



Duvernell, D. D., J. F. Schaefer, D. C. Hancks, ]. A. Fonoti and A. 
M. Ravanelli. 2007. Hybridization and reproductive isolation 
among syntopic populations of the topminnows Fundulus 
notatus and F. olivaceus. J. Evol. Biol. 20:152-164. 

Endler, J. 1993. The color of light in forests and its implications. 
Ecol. Monogr. 63:1-27. 

Fraser, D. E, J. F. Gilliam, and T. Yip Hoi. 1995. Predation as an 
agent of population fragmentation in a tropical watershed. 
Ecology 76:1461-1472. 

Fuller, R. C, and ). Travis. 2004. Genetics, lighting environment, 
and heritable responses to lighting environment affect male 
color morph expression in bluefin killifish, Lucania goodei. 
Evolution 58:1086-1098. 

Fuller, R. C, K. L. Carleton, J. M. Fadool, T. C. Spady, and ]. 
Travis. 2005. Genetic and environmental variation in the visual 
properties of bluefin killifish, Lucania goodei. J. Evol. Biol. 
18:516-523. 

Fuller, R. C, K. E. McGhee, and M. Schrader. 2007. Speciation in 
killifish and the role of salt tolerance. J. Evol. Biol. 
20:1962-1975. 

Gido, K. B., J. Schaefer, and J. A. Falke. 2009. Convergence of fish 
communities from the littoral zone of reservoirs. Freshw. Biol. 
54:1163-1177. 

Gordon, S. P., A. Lopez-Sepulcre, and D. N. Reznick. 2011. 
Predation-associated differences in sex-linkage of wild guppy 
coloration. Evolution 66:912-918. 

Gutierrez, M. 2010. Reproductive behaviors of male and female 
blackspotted topminnows, Fundulus olivaceus. MS thesis, The 
University of Southern Mississippi. Hattiesburg, MS, p. 38. 

Haskins, C. P., E. F. Haskins, J. J. A. McLaughing, and R. E. Hewit. 
1961. Polymorphisms and population structure in Lebistes 
reticulates, an ecological study. Pp. 320-395 in W. F. Blair, ed. 
Vertebrate speciation. Univ. of Texas Press, Austin, TX. 

Hendry, A. P., M. L. Kelly, M. T. Kinnison, and D. N. Reznick. 
2006. Parallel evolution of the sexes? Effects of predation and 
habitat features on the size and shape of wild guppies. J. Evol. 
Biol. 19:741-754. 

Houde, A. E., and }. A. Endler. 1990. Correlated evolution of 
female mating preferences and male color patterns in the 
guppy Poecilia reticulata. Science 248:1405-1408. 

Howell, W. M., and A. Black. 1981. Karyotypes in populations of 
the cyprinodontid fishes of the Fundulus notatus species 
complex: a geographic analysis. Bull. Ala. Mus. Nat. Hist. 
6:19-30. 

Ingleby, F. C, ]. Hunt, and D. J. Hosken. 2010. The role of 

genotype-by-environment interactions in sexual selection. 

J. Evol. Biol. 23:2031-2045. 
Iwasa, Y., and A. Pomiankowski. 1999. Good parent and good 

genes models of handicap evolution. J. Theor. Biol. 

200:97-109. 

Iwasa, Y, A. Pomiankowski, and S. Nee. 1991. The Evolution of 
costly mate preferences II. The "Handicap" principle. 
Evolution 45:1431-1442. 

Knight, G. L., and K. B. Gido. 2005. Habitat use and susceptibility 
to predation of four prairie stream fishes: implications for 



1380 



© 2012 The Authors. Published by Blackwell Publishing Ltd. 



J. F. Schaefer et al. 

conservation of the endangered Topeka shiner. Copeia 
2005:38-47. 

Krebs, J. M., and S. S. Bell. 2011. Risk of predation reflects 

variation in the reproductive strategy of a dominant forage fish 
in mangrove tidal tributaries. Oecologia. 168:737-748. 

Langerhans, R. B. 2008. Predictability of phenotypic 

differentiation across flow regimes in fishes. Integr. Comp. 
Biol. 48:750-768. 

Langerhans, R. B., and D. Reznick. 2010. Ecology and evolution 
of swimming performance in fishes: predicting evolution with 
biomechanics. Pp. 200-248 in P. Domenici and B. G. Kapoor, 
eds. Fish locomotion: an eco-ethological perspective. Science 
Publishers, Enfield, New Hampshire. 

Langerhans, R. B., C. A. Layman, A. Shokrollahi, and T. J. 
DeWitt. 2004. Predator-driven phenotypic diversification in 
Gambusia affinis. Evolution 58:2305-2318. 

Lee, D. S., C. R. Gilbert, C. H. Hocutt, R. E. Jenkins, D. E. 
McAllister, and J. R. Stauffer. 1980. Atlas of North American 
freshwater fishes. North Carolina Biological Survey. Raleigh, 
NC. 

Lindholm, A., and F. Breden. 2002. Sex chromosomes and sexual 
selection in poeciliid fishes. Am. Nat. 160:S214-S224. 

Matthews, W. J. 1998. Patterns in freshwater fish ecology. 
Chapman and Hall, New York, NY. 

Matthews, W. J., K. B. Gido, G. P. Garrett, F. P. Gelwick, J. G. 
Stewart, and J. F. Schaefer. 2006. Modular experimental 
riffle-pool stream system. Trans. Am. Fish. Soc. 
135:1559-1566. 

Millar, N. P., D. Reznick, M. T. Kinnison, and A. P. Hendry. 2006. 

Disentangling the selective factors that act on male colour in 

wild guppies. Oikos 113:1-12. 
Moller, A. P., and R. V. Alatalo. 1999. Good-genes effects in 

sexual selection. Proc. R. Soc. Lond. B 266:85-91. 
Parker, T. H., R. Knapp, and J. A. Rosenfield. 2011. Social 

mediation of sexually selected ornamentation and steroid 

hormone levels in male junglefowl. Anim. Behav. 64:291-298. 
R Development Core Team. 2010. R: a language and environment 

for statistical computing. R Foundation for Statistical 

Computing, Vienna, Austria. 
Reznick, D., and J. A. Endler. 1982. The impact of predation on 

life history evolution in Trinidadian guppies (Poecilia 

reticulata). Evolution 36:160-177. 
Reznick, D., M. J. Butler, and H. Rodd. 2001. Life-history 

evolution in guppies. VII. The comparative ecology of high- 

and low-predation environments. Am. Nat. 157:126-140. 
Rosenthal, G. G., T. Y. Flores Martinez, F. J. Garcia de Leon, and 

M. J. Ryan. 2001. Shared preferences by predators and females 

for male ornaments in swordtails. Am. Nat. 158:146-154. 
Ross, S. T. 2000. Inland fishes of Mississippi. Univ. Press of 

Mississippi, Oxford, MS. 
Schaefer, J., B. R. Kreiser, C. Champagne, P. M. Mickle, and D. D. 

Duvernell. 2009. Patterns of co-existence and hybridisation 

between narrowly endemic (Fundulus euryzonus) and broadly 

distributed (F. olivaceus) topminnows in a riverine contact 

zone. Ecol. Freshw. Fish. 18:360-368. 



Evolution of Sexual Dimorphism in Topminnows 

Schaefer, J. F., D. D. Duvernell, and B. R. Kreiser. 2011. Ecological 
and genetic assessment of spatial structure among replicate 
contact zones between two topminnow species. Evol. Ecol. 
25:1145-1161. 

Setzer, P. Y. 1970. An analysis of a natural hybrid swarm by means 
of chromosome morphology. Trans. Am. Fish. Soc. 
99:139-146. 

Sharma, M. D., T. Tregenza, and D. J. Hosken. 2010. Female mate 

preferences in Drosophila simulans: evolution and costs. J. Evol. 

Biol. 23:1672-1679. 
Stevens, M., and S. Merilaita. 2009. Defining disruptive 

coloration and distinguishing its functions. Philos. Trans. R. 

Soc. B 364:481-488. 
Suttkus, R. D., and R. C. Cashner. 1981. A new species of 

cyprinodontid fish, genus Fundulus (Zygonectes), from Lake 

Pontchartrain tributaries in Louisiana and Mississippi. Bull. 

Ala. Mus. Nat. Hist. 6:1-17. 
Sutton, S. G., T. P. Bult, and R. L. Haedrich. 2000. Relationships 

among fat weight, body weight, water weight, and condition 

factors in wild Atlantic salmon parr. Trans. Am. Fish. Soc. 

129:527-538. 

Thomerson, J. E. 1966. A comparative biosystematic 
study of Fundulus notatus and Fundulus olivaceus 
(Pisces: Cyprinodontidae). Tulane Stud. Zool. 13:29— 
47. 

Thomerson, J. E. 1967. Hybrids between the cyprinodontid 
fishes, Fundulus notatus and Fundulus olivaceus in Southern 
Illinois. Trans. Illi. St. Acad. Sci. 60:375-379. 

Thomerson, J. E., and D. P. Woolridge. 1970. Food habits of 
allotopic and syntopic populations of the topminnows 
Fundulus olivaceus and Fundulus notatus. Am. Mid. Nat. 
84:573-576. 

Vannote, R. L., G. W. Minshall, K. W. Cummins, J. R. Sedell, and 

C. E. Cushing. 1980. The river continuum concept. Can. J. 

Fish. Aq. Sci. 37:130-137. 
Vergara, P., J. Martinez-Padilla, F. Mougeot, F. Leckie, and S. M. 

Redpath. 2011. Environmental heterogeneity influences the 

reliability of secondary sexual traits as condition indicators. 

J. Evol. Biol. 35:20-28. 
Vigueira, P., J. F. Schaefer, D. D. Duvernell, and B. R. Kreiser. 

2008. Tests of reproductive isolation among species in the 

Fundulus notatus (Cyprinodontiformes: Fundulidae) species 

complex. Evol. Ecol. 22:55-70. 
Walsh, M. R., and D. N. Reznick. 2009. Phenotypic diversification 

across an environmental gradient: a role for predators and 

resource availability on the evolution of life histories. 

Evolution 63:3201-3213. 
Walsh, M. R., and D. N. Reznick. 2010. Influence of the indirect 

effects of guppies on life-history evolution in Rivulus hartii. 

Evolution 64:1583-1593. 
West-Eberhard, M. J. 2003. Developmental plasticity and 

evolution. Oxford Univ. Press, New York. 
Whitehead, A. 2010. The evolutionary radiation of diverse 

osmotolerant physiologies in killifish (Fundulus sp.). Evolution 

64:2070-2085. 



©2012 The Authors. Published by Blackwell Publishing Ltd. 



1381 



